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SUMMARY 


A broad base of thermocouple and phase-change-paint data was assembled and 
correlated to the nominal-design 14414.1 and proposed STS-1 (first flight of the 
space transportation system) entry trajectories. Averaged data from phase-change- 
paint tests compared favorably with thermocouple data for predicting heating rates. 
Laminar and turbulent radiation-equilibrium heating rates were computed on the lower 
surface of the Shuttle orbiter for both trajectories, and the lower-surface center- 
line results were compared both with aerodynamic-heating design data and with flight 
values from the STS-1 and STS-2 trajectories. The peak laminar-heating values from 
the aerodynamic-heating design-data book were generally 40 to 60 percent higher than 
the laminar estimates of this study, except at the 55-percent location of maximum 
span where the design-data-book values were less than 10 percent higher. 

Estimates of both laminar and turbulent heating rates compared favorably with 
flight data. Peak turbulent heating rates were estimated to occur 200 to 400 sec 
later in the trajectory than the laminar peaks and could be more than twice the 
magnitude of the design-data values. Therefore, a precise determination of flight 
boundary-layer transition becomes a critical exercise in thermal-protection-system 
evaluation. The cause of the increased heating in the vicinity of 55 percent of the 
span could not be resolved with the simplified approach employed in this study. 
However, the phase-change-paint data were found to be a very effective means of 
mapping the region of influence. 


INTRODUCTION 

An aerothermal analysis of the heating on the lower surface of the Shuttle 
orbiter during entry was initiated at the Langley Research Center as direct support 
to a study sponsored by the Johnson Space Center of the strength and integrity of the 
thermal protection system (TPS) of the Space Shuttle orbiter. The analysis reported 
herein was intended to be an independent assessment of the radiation-equilibrium 
heating rates and surface temperatures expected on the lower surface of the orbiter 
for the nominal-design trajectory. In order to design the primary TPS, these results 
would then be compared with the data contained in "Space Shuttle Orbiter Entry 
Aerodynamic Heating Data Book" by J. W. Haney and C. T. Petrilla (SD73-SH-01 84, 

C Revision, Book I, October 1978) which were established by the Rockwell Inter- 
national Corporation. The existing base of experimental tunnel heat-transfer data 
was extensively surveyed to collect the most appropriate data to form a sound 
empirical approach to heating predictions. The resulting correlations for laminar 
heating rates are a vital part of the approach and were formulated on the basis of 
combined thermocouple and phase-change-paint data. 

The estimates of local inviscid conditions and flight stagnation-point heating 
rate required to apply the correlations were generated by straightforward and well- 
known techniques. The method of Fay and Riddell (ref. 1) was used to determine 
stagnation-point heating rates. Local inviscid-f low conditions were based on a 
modification of a tangent-cone approximation technique presented in reference 2. 



The turbulent heating rates presented in this paper were generated by laminar- 
to-turbulent heat-transfer ratios based on well-known approximations to classical 
boundary-layer equations . The same approach was also subsequently applied in the 
Shuttle Flight Certification Program. During flight certification, the analysis was 
extended to include portions of the orbiter other than the lower surface; however, 
the results presented herein contain only lower-surface data. Neither the empirical 
approach developed and employed for these studies nor the predicted results have 
appeared in a referenceable document. 


SYMBOLS 

In addition to a technique description and a discussion of calculated results, 
data from the lower-surface center line are compared with flight data from the first 
and second flights of the space transportation systems (STS-1 and STS-2) to assess 
both the analytical approach and predicted results . 

b wing span 

C pressure coefficient; specific heat 

P 

c local axial chord length 

2 

g gravitational constant, 32.152 ft/sec 

h altitude 

\ body length 

M Mach number 

effective Mach number normal to shock 

Prandtl number 
Stanton number 
static pressure 
heating rate per unit area 
convective heating rate per unit area 

Fay-Riddell radiation-equilibrium heating rate per unit area at 
stagnation point of 1-ft-radius sphere in flight 

R local Reynolds number in conical flow 

O / X 

R* reference- temperature Reynolds number 

R^ ^ free-stream Reynolds number based on orbiter body length 

T temperature 

t time 
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u 

film heat-transfer coefficient 


^ref 

film heat-transfer coefficient at stagnation point of 1-ft-radius sphe 
flight (or its scale equivalent) 

V 

velocity 


X 

axial distance aft of leading "edge (or nose) 


y 

lateral distance spanwise from fuselage center 

line 

a 

angle of attack 


Y 

ratio of specific heats for air 


e 

surface emissivity 


\ 

recovery factor 


0 

s 

shock angle 


M. 

viscosity coefficient 


P 

density 


a 

Stef an-Boltzmann constant for radiant-heating 

calculations 

Subscripts; 


aw 

adiabatic wall 


c 

conical-flow conditions 


lam 

laminar 


o 

stagnation or total conditions 


pc 

phase-change paint 


ref 

radiation-equilibrium conditions at stagnation 
sphere in flight (or its scale equivalent) 

point of 1-ft-radius 

s 

conditions behind oblique shock 


turb 

turbulent 


w 

wall 


CO 

free stream 


Superscript: 


* 

reference-temperature conditions 
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Abbreviations : 


AEDC Arnold Engineering Development Center 

STS space transportation system 

TPS thermal protection system 


ANALYTICAL APPROACH 

The most desirable and comprehensive approach to the aerothermal-analysis task 
would be a coupled inviscid and viscous theoretical solution of either the complete 
vehicle or the windward-side flow field at the entry flight conditions representing 
peak heating coupled with appropriate validating experimental results. This is 
particularly true if the basic solutions could be systematically perturbed and 
validated by experimental data in order to identify and evaluate parameters which 
most significantly influence the orbiter's aerothermal environment. 

Areas of concern for the lower surface include real-gas effects (equilibrium and 
finite-rate chemistries), geometric influences such as abrupt changes in leading-edge 
sweep and control surfaces, shock interactions, shear layers downstream of shock 
interactions, boundary-layer transition, and excursions in flight profile. The sides 
and upper surfaces may have additional complexities resulting from separation, inter- 
acting or interfering flows, and possible vortex scrubbing. After surveying the 
available numerical techniques, however, it was quickly determined that a compre- 
hensive theoretical approach was not practical for this effort. This decision was 
reached because a good coupled inviscid and viscous solution for complex configura- 
tions like the orbiter at its intended flight conditions was not state of the art; 
also, the probability of successfully applying such codes within the prohibitively 
short time frame of either the TPS strength-and-integrity study or the first flight 
certification would have been very tenuous and could have jeopardized the possibility 
of any real assessment of the heating. 

Therefore, an alternative approach based on the existing Shuttle aerothermo- 
dynamic data base was adopted as the best approach for quickly meeting both the 
TPS strength-and-integrity study objectives and subsequent first-flight-certification 
objectives. A survey of available data was conducted to collect the most pertinent 
heat-transfer data. Contacts were also made with several past experimental analysts 
who were active during the development phases of the Shuttle to secure additional 
insight and data. The most consistent and coitprehensive heating data generated in 
any one facility seemed to be the data obtained at a Mach number of 8 in the AEDC 
von Karman Gas Dynamics Facility Tunnel B (referred to herein as AEDC Tunnel B) . 

Three series of heating tests (refs. 3, 4, and 5) were conducted on 0.0175-scale 
models at angles of attack and Reynolds numbers based on length which could be used 
to approximate key portions of the nominal- design 14414.1 entry trajectory illus- 
trated in figure 1(a) as well as the proposed STS-1 trajectory illustrated in fig- 
ure 1(b). These data do not contain real-gas effects. Therefore, the predicted 
heating rates could be significantly higher than those measured in flight in the 
early portions of the trajectories where considerable nonequilibrium flows are known 
to exist; however, as the flows tend to reach equilibrium later in the trajectories, 
the real-gas effects should tend to become negligible. (See ref. 6. ) 

The initial 1450-sec interval of each trajectory covered all the major heating 
conditions and was the basic flight interval used throughout this study. The peak 
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heating condition usually occurs between 400 and 1000 sec as illustrated by the 

curves which represent the stagnation-point radiation-equilibrium heating rate 
on a 1-ft-radius sphere. (See fig. 1.) The two trajectories of interest are defined 
in figures 2 and 3. Four sets of experimental test conditions were used to approxi- 
mate the flight Reynolds numbers and angles of attack of the nominal-design tra- 
jectory as shown in figures 2(a) and 3(a), respectively. Five corresponding stepwise 
approximations for the proposed STS-1 entry trajectory are presented in figures 2(b) 
and 3(b). No interpolations of the experimental data were performed during this 
analysis, and the step functions illustrated in figures 2 and 3 cause slight discon- 
tinuities that are sometimes evident in subsequent heating histories. 

Four spanwise locations were selected to represent the key features of the lower 
surface based on a visual analysis of all experimental data. The locations include 
the lower center line and the 40- , 55-, and 70-percent locations of maximum span. No 
thermocouples were located at the 55-percent span location; however, the phase-change- 
paint data consistently indicated a narrow region of increased heating there. The 
thermocouple data obtained along the 50- and 60-percent spans do not indicate 
increased heating. 


Thermocouple Data 

The thermocouple pattern and test conditions for the thermocouple data employed 
are illustrated in figure 4, There were an additional 23 thermocouples on the for- 
ward portion of the lower center line which are not shown. The basic form of the 
data employed is the film heat-transfer coefficient U nondimensionalized by the 
coefficient for stagnation-point radiation-equilibrium heating to a scale equivalent 
to a 1-ft-radius sphere in flight. The equivalent wind-tunnel sphere radius was 
0.0175 ft in references 3, 4, and 5. 


Phase-Change-Paint Data 

A fairly comprehensive description of the phase-change-paint test technique and 
the assumptions required to obtain aerodynamic-heating estimates are contained in 
reference 7. Basically, a model constructed of a dark-colored insulating material 
having well-defined thermal properties which are well-known is coated with a 
temperature-sensitive paint prior to injection into an airstream. The paint employed 
is a light-colored opaque solid which becomes a colorless liquid at a specific tem- 
perature. Paints which change phase at almost any given temperature are available. 
When the model is injected into the test airstream, the paint in the areas of high 
heating melts first and exposes the darker model surface with the peripheries of 
these areas being denoted as melt lines. These areas, which enlarge with time, are 
photographed with a motion-picture camera operating at a known preset rate. These 
photographs give a precise time history of the distinct melt lines which, coupled 
with a knowledge of the model thermal properties, presents sufficient information to 
conduct a one-dimensional transient heat-conduction analysis and to extract surface 
heat- transfer coefficients. 

Tracings of the phase-change-paint melt lines and of the heat- transfer- 
coefficient ratios U/U^gf that they represent, approximating the earliest portions 
of both the design (14414.1) and STS-1 trajectories, are illustrated in figure 5 for 
an angle of attack of 40° and a Reynolds number of 1 x 10^. Figure 5(a) was gener- 
ated from the data presented in reference 4. A key objective of that study was to 
investigate the heating of a series of deflected elevens; however, only those data 
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representing undeflected elevens were employed here. Figures 5(b) to 5(j) were 
generated from data contained in reference 5. A major objective of that study was 
to investigate blockage and the possible influence of local cross flows that might be 
established by cutting away portions (notches) of one side of the model as evident in 
figures 5(h) to 5(j). The data from that study compared favorably with the complete 
orbiter data to the extent that no apparent effects due to the notches were observed. 
Figure 6 (taken from ref. 4) contains the only set of phase-change-paint data 
applicable to the plateau denoted by a = 35° and « 1 .9 x 10°, which was used 

to define the second time interval of the nominal trajectory where 
assumed constant. The interval denoted by a » 30° and R « 3.75 x l 0® is 
represented by figure 7. The data shown in figures 7(a) to°°')^d) were taken from 
reference 4, and the data shown in figures 7(e) to 7(q) were taken from reference 5. 
Phase-change-paint data applicable to the last step function of both trajectories 
(a « 30° and « 7 x 10^) are represented in figures 8(a) and 8(b), which were 

derived from reference 4, and in figures 8(c) to 8(j), which were generated from data 
contained in reference 5. 


Averaged Thermocouple and Phase-Change-Paint Data 

Values of for the 0-, 40-, 55-, and 70-percent locations of maximum 

span were extracted from the data sets and plotted as functions of nondimens ional 
axial location x/c. The data for constant angle of attack and Reynolds number that 
were used to represent both the first 600 sec of the nominal-design 14414.1 tra- 
jectory and the first 850 sec of the STS-1 trajectory are illustrated in figure 9. 

The circular symbols always represent the thermocouple data so that an easy com- 
parison can be made with the phase-change-paint data. The symbol key shown in the 
upper right corner of figure 9 is consistent with those provided at the end of the 
test-condition legends in figure 5 so that individual phase-change-paint tracings 
can be easily keyed to the data plots. The curves represent the author's "simple 
fairing" of the individual data sets. Data from these curves, rather than the indi- 
vidual data points, were used in subsequent heating calculations along the entry 
trajectory for direct comparison with the aerodynamic-heating data book. A more rig- 
orous and systematic fairing of the data could be based on laminar and/or turbulent 
decays in heating with distance. 

The experimental data to approximate the trajectory parameters between 600 and 
850 sec used an angle of attack of 35° and a Reynolds number based on length of 
approximately 2 x 1 0® as illustrated in figure 10. The scarcity of data available 
for a = 35° makes the curve fairing at this condition questionable. The experi- 
mental data sets used to approximate between 850 and 1 050 sec of the nominal-design 
14414.1 trajectory are illustrated in figure 11. The data of figure 12 were used 
over the interval between 1050 and 1450 sec of the nominal-design 14414.1 trajectory 
and also between 1250 and 1450 sec for the STS-1 entry trajectory. Both of these 
data sets were obtained at a = 30° but at separate Reynolds numbers. The data 
obtained at a Reynolds number of 4 x i and shown in figure 11(a) appear laminar 
along the lower center line. However, a comparison of these data with the low 
Reynolds number data of figure 9 indicates probable boundary-layer transition at the 
55-percent span location and possibly at the 40- and 70-percent span locations as 
well. The highest Reynolds number data (R^ ^ » 8 x i o^ shown in fig. 12) exhibit 
boundary-layer transition at all spanwise locations. Transition seems to occur about 
55 percent of the axial distance down the lower center line, about 40 to 45 percent 
of the chord at 2y/b = 0.40, about 15 percent of the chord at 2y/b = 0.55, and 
about 20 to 25 percent of the chord at 2y/b = 0.70. 
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Figures 13, 14, and 15 contain the additional averaged experimental data which 
were combined with the data in figures 9 and 12 to approximate the proposed STS-1 
trajectory. The heat-transfer-coefficient ratios presented in figure 13 (a = 40° 
and ^ « 2 x 10®), representing the interval between 850 and 1000 sec of the pro- 

posed sfe-l entry, appear to be essentially laminar at all four spanwise locations. 
The data in figure 14 representing the next interval between 1000 and 1150 sec are 
also at a = 40 ° but at approximately twice the Reynolds number. There is an indi- 
cation of possible transition at the 0-, 40-, and 55-percent spanwise locations, 
whereas the most outboard location has a more laminar characteristic longitudinal 
distribution. The data representing the STS-1 interval between 1150 and 1250 sec 
(a = 35° and ^ « 8 x 10®) clearly exhibit transition at all four spanwise loca- 
tions. Transition seems to occur at about 60 percent of the local chord for the 
center line and at about 30 percent of the local chord at the other three spanwise 
locations . 


METHODOLOGY FOR EXTRAPOLATING WIND-TUNNEL DATA TO FLIGHT 

The heat-transfer-coefficient ratios measured at a Mach number of 8 in the AEDC 
Tunnel B at flight angles of attack and flight free-stream Reynolds numbers based on 
vehicle length were assumed to represent adequately the ratios expected in those 
portions of the flight trajectory not dominated by real-gas effects. No corrections 
for Mach number and wall-temperature effects were made. A simple procedure for 
calculating flight radiation-equilibrium heating rates and surface temperatures was 
derived in the following manner: 

If the convective heating is defined by 


q = U(T - T ) 
^conv aw w 


(1 ) 


then multiplying both sides of equation (1) by ^ref^'^ref 
gives 


and collecting terms 


q 

conv 


q 


ref 


U 


T - T 
aw w 


u 


ref 


ref 


ref 


( 2 ) 


The subscript ref refers to stagnation-point heating on a 1-ft-radius sphere or a 
scaled equivalent. The relationship for convective heating to the nose of the sphere 
can be defined as 


'^ref ^ref^'^o "^ref^ 


Similarly, the reference-condition relation for radiant heating is 


q 


ref 


aeT 


4 

ref 


(3) 


(4) 
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Thus, is defined as a radiation-equilibrium heating rate. Therefore, solving 

equation (4) for substituting Tj.^^ into equation (3), and collecting terms 

yields 


U 

ref 



(5) 


Substituting this relation for into equation (2) and simplifying gives 


q 

conv 


U 


ref U 


ref 



(6) 


Setting the right-hand side of equation (6) equal to the radiant-heating relation 
gives 



U 


ref U 


ref 


T - T 
aw w 



(7) 


which can be evaluated for radiation-equilibrium wall temperatures at flight condi- 
tions by a simple iterative procedure. The emissivity e was held constant at 0.9 
throughout. 


Methodology for Calculating Turbulent Heating 

The procedure used to calculate turbulent heating was based on conical-flow edge 
conditions. The conical-flow assumption was adopted on the basis of previous work at 
the Johnson Space Center as illustrated in figure 16 where experimental pressure data 
along the lower-surface center line are compared with wedge and conical-flow esti- 
mates. The conical-edge conditions were used in conjunction with empirical reference- 
temperature techniques to calculate a turbulent-to-laminar heat-transfer-coefficient 
ratio ^turb/^lam* turbulent flight heating was then determined by simply apply- 

ing this ratio to the laminar flight values based on the wind-tunnel data. 

The following approach was established to approximate conical-flow conditions at 
the edge of the boundary layer. A conical shock angle was estimated from the simple 
relation 


0 

s 


1 .035 
a 


( 8 ) 
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which is an empirical fit of data from reference 8 and gives reasonable results for 
Mach numbers above 10. The following relations were also taken from reference 8 to 
calculate conditions immediately behind the conical shock angle 9^; 



2vM^ sin^S 
' oo s 

T+ 


(y - 1) 


(9) 



2 ( 1 ^ - 1 ) 

M^[(y + 1)5 + (y - 1)] 


( 10 ) 


^ = r (y-'i)C+(y + 'i) 

^[(y + 1)5 + (y - 1). 


( 11 ) 



( 12 ) 



S 


M^[(y + 1)5 + (y - 1)] - 2(5^ - 1) 

5[(y - 1)5 + (y + Dl 


(13) 


The conical pressure at the boundary-layer edge was approximated by 




where is an effective Mach number normal to the shock and is defined by 


2 . 2 
yM sin a 



M = (0.87M - 0.554) sin a - 0.53 

ns “ 


( 15 ) 
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This expression coupled with equation (14) has been assessed and shown to be very 
accurate in an appendix to reference 9. The conical compression between shock and 
boundary-layer edge was calculated from equations (9) and (14): 




(16) 


The following isentropic relations were employed to define the edge conditions 
further: 


Jzl 




(17) 


(18) 


(19) 


Viscosity was determined from the Sutherland equation 


2.329 X 10 
c 

^c ■ T + 215.933 
c 


and local Mach number was determined from 


c 


P 

c c 


A set of modified Eckert reference temperatures were calculated by 


* 

— = 1 .0 + 0.0393n 
T r c 

c 


+ 0 



( 20 ) 


(21 ) 


( 22 ) 
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where is the wall temperatur e and t)^ is the recovery factor which was taken to 

be ^Npj, for laminar flow and for turbulent flow. The Prandtl number Np^ 

was held constant at 0.71 throughout this study. A reference Reynolds number was 
calculated from 


* 

R = 

X 



(23) 


where 


R 

c,x 


p V X 
^c c 

^^c 


(24) 


A laminar Stanton nxamber was approximated from 


N 


St, lam 


0.332\j3 



(25) 


and a turbulent Stanton number was approximated from 


N 


0.222(1 .15) 


St, turb 


(109, 0 4) 


2.584 


(26) 


which is based on a Prandtl-Schlichting relation for average skin-friction coeffi- 
cient. Since the heat-transfer coefficient can be defined in terms of Stanton number 
as 


u = gvpCpNg^ 

the turbulent heat-transfer coefficient becomes simply 


(27) 


N 


U 


= u 


turb lam N 


St, turb 
St, lam 


The turbulent convective heating can be defined as 


q = , (T - T ) 

turb turb aw w 


(28) 


( 29 ) 



and set equal to the radiant heating to yield 


aeT^ = - T ) (30) 

w turb aw w 

Equations (22) to (30) can be iterated for the turbulent radiation-equilibrium wall 
temperature. One or two iterations are all that are required if the values of 
laminar radiation-equilibrium temperatures are used as a first guess. 


RESULTS AND DISCUSSION 

Comparisons for Nominal-Design 14414.1 Entry Trajectory 

Laminar heating .- The experimental values of U/U^^^ were assumed laminar and 
were utilized to calculate the laminar radiation-equilibrium heating rates and sur- 
face temperatures at four spanwise locations over the nominal-design 14414.1 entry 
trajectory. Figure 17 presents a summary comparison of the peak laminar heating 
rates from the aerodynamic-heating design-data book with those determined in this 

study which occur at R » 1 x 10®. The present estimates are considerably lower 

1 

than those from the design-data book, approximately 60 percent at the center line and 
the outboard span location (2y/b = 0.70). In the region beginning at 2y/b = 0.40 
and ending at approximately 2y/b = 0.55, the conservatism in peak heating exhibited 
by the design-data book is diminished. This diminished conservatism is most evident 
at the 55-percent span location where less than a 10-percent margin is exhibited 
beginning approximately halfway down the wing chord (x/c = 0.5). However, it should 
be noted that the design data are still conservative everywhere when compared with 
the experimental laminar-heating predictions presented here. 

Turbulent heating.- An estimate of turbulent heating rates was also made for the 
nominal-design 14414.1 trajectory. This analysis did not depend on any estimate of 
onset of transition; instead, a turbulent-heating calculation was made by scaling the 
laminar-heating values of the wind-tunnel data for throughout the trajectory 

to define an upper boundary for the heating. In this way the effect of any subse- 
quent estimate of transition on lower-surface heating could be evaluated immediately. 

Comparison of design and current estimates .- Figure 18 presents the present- 
study estimates of laminar and turbulent heating-rate histories as well as those from 
the aerodynamic-heating design-data book at five axial locations (x/c = 0.10, 0.25, 
0.50, 0.75, and 0.90) for the four locations of maximimi span (2y/b =0.00, 0.40, 

0.55, and 0.70). The 10-percent axial-location (x/c = 0.10) data are illustrated in 
figure 18(a). The solid curves are direct traces of the design-data-book histories. 
The data-book four-digit identification numbers (body points), as well as scales and 
units of heating rates, were retained for easy reference. The design-data-book four- 
digit curve letters (ABCD) can also be used to identify quickly the general surface 
locations for the data presented herein. When A = 1 and D = 0 (i.e., 1BC0), the 

curve is identified with the lower-surface fuselage center line. For the center-line 
curves, BC defines the nondimensional axial location x/\. When A = 2, the curve 
corresponds to a location on the lower surface of the wing. For these lower-wing 
curves, B identifies the orbiter half-span location (i.e,, B = 2y/b x 10) and CD 
identifies the percentage of local wing chord (i.e., CD = x/c x 100). The design- 
data book does not contain curves at the 55-percent spanwise location; therefore, 
the curves at the 50-percent (2510) and the 60-percent (2610) span locations were 



included for comparison. The jagged appearance of the present calculations is 
primarily a result of the stepwise approximation of the design trajectory. 

The heating estimates from the aerodynamic-heating design-data book are suffi- 
cient to accommodate turbulent heating for the forward locations (x/c = 0.10) on 
the 40- and 70-percent span locations. At the center line and the 55-percent span 
locations, however, the current estimate of peak turbulent heating is approximately 
30 percent higher than the design data. Greater disparity (turbulent to laminar and 
turbulent to design-data book) occurs at increasing axial distances from the nose and 
leading edges as illustrated in figures 18(b) and 18(e). When the design-data book 
did not contain body points in close proximity to the calculated values, the curves 
on either side were reproduced and presented for comparison. For example, at 
25-percent chord for the 55-percent span location, the comparison included axial 
points before (xx20) and after (xx30) as well as spans inboard (x5xx) and outboard 
(x6xx) of the point of interest. 

The peak turbulent heating occurs at a later time in the trajectory than the 
laminar peaks which is in the correct direction to coincide with the changes in edge 
and reference Reynolds numbers. The reasonably good agreement in slopes and magni- 
tudes between the turbulent heating-rate histories from this study and the turbulent- 
heating portions of the design-data-book curves illustrated in figure 18 is gratify- 
ing; however, the simplified approach employed here may not be rigorous enough to 
exclude fortuity. However, a partial validity check of the turblent-heating method- 
ology was conducted by using the method of reference 10 at t = 480 sec to calculate 
the equivalent equilibrium airflow conditions corresponding to perfect-gas parameters 
defined by equations (8) to (22). Even though large variations occurred in some 
parameters, like static temperature and local Mach number, the ratio of turbulent to 
laminar Stanton numbers increased less than 10 percent. Therefore, the analysis 
appears adequate to infer that in the event of possible transition at 1000 sec (or 
earlier) into the entry, a more comprehensive and rigorous effort may be required to 
validate the TPS. 


Heating Estimates for Proposed STS-1 Entry Trajectory 

During the Shuttle Flight Certification Program the same methodology was 
employed to generate rapid independent estimates of the aerothermal environment 
associated with the proposed STS-1 entry trajectory depicted graphically in fig- 
ure 1(b). Although only results for the lower surface are presented here, many other 
regions on the orbiter were examined during the certification prior to first flight. 
The upper wing, fuselage sides and top, canopy region, orbiter-maneuvering-system 
(QMS) pods, vertical tail, and local surface discontinuities such as the reaction- 
control-system (RCS) nozzles and payload bay joints and hinges were examined by this 
technique. The results presented here are limited to the same regions previously 
discussed for the TPS strength-and-integrity study: namely, the lower-surface center 

line along with the 40-, 55-, and 70-percent locations of maximum span. The five 
sets of averaged experimental heat-transfer-coefficient ratios employed to approxi- 
mate the proposed STS-1 trajectory are presented in figures 9 and 12 to 15. 

Since the proposed STS-1 trajectory was expected to produce a somewhat more 
benign aerothermal environment than the nominal-design 14414.1 trajectory, the 
laminar- and turbulent-heating estimates for both trajectories are presented together 
for easy comparison in figure 19. A general reduction in peak laminar heating during 
the STS-1 trajectory for longitudinal locations less than 75 percent of the local 
chord at all spanwise locations is readily apparent in figures 19(a) to 19(d). The 
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differences in heating rates and any associated aerothermal advantages during the 
STS-1 trajectory become less clear or pronounced toward the rear of the orbiter, 
probably because of the transitional and/or turbulent wind-tunnel data that had to 
be employed. 

The estimated turbulent boundary peaks for the two trajectories tend to be more 
nearly equal than the laminar peaks. The major difference seems to be a later occur- 
rence of the turbulent peaks during the STS-1 trajectory than during the nominal- 
entry 14414.1 trajectory. However, based on a comparison of some key parameters, for 
example, the Reynolds number curves presented in figure 2, it is likely that transi- 
tion is also delayed approximately the same time interval. In any event, the values 
from the aerodynamic-heating design-data book for peak laminar heating are adequate 
to accommodate turbulent heating if it should occur later than 1000 sec into the 
nominal-design 14414.1 entry trajectory, and into the STS-1 entry as well, as long 
as transition (fully turbulent flow) does not occur prior to 1200 sec into the entry 
trajectory. 


Data Comparisons of Flight Heating Rates 

Flight heating rates along the lower-surface center line of the orbiter measured 
during the STS-1 arid STS-2 entries were presented for a series of flight times and 
conditions in reference 11. Key trajectory parameters at eight times during the 
STS-1 entry and seven times during the STS-2 entry are compared with the nominal- 
design 14414.1 and the proposed STS-1 trajectories in figure 20. The two flight 
histories are fairly well represented by the proposed STS-1 trajectory. However, the 
Reynolds number, which is a key heating parameter, tends to be slightly higher than 
that proposed with the result that between 1100 and 1300 sec into the entry, it more 
nearly matches the Reynolds number predicted for the nominal-design 14414.1 entry 
trajectory. 

The heating rates along the lower-surface center line are compared in figure 21 
at the 10-, 25-, 50-, 75-, and 90-percent locations of the fuselage length. The 
laminar and turbulent heating rates estimated for the proposed STS-1 trajectory, as 
well as the values from the aerodynamic-heating design-data book, are presented for 
each longitudinal location. Except for the earliest STS-2 data point corresponding 
to 500 sec on the forward portion of the orbiter (x/c < 0.50) where nonequilibrium 
flow might be expected, the estimates appear to compare very favorably with the 
flight measurements. Fortunately, when disparities are noted, the estimated values 
tend to be conservative. 

The heat-transfer-coefficient ratios from reference 3 were cross-plotted as 
functions of Reynolds number and angle of attack to determine wind-tunnel transition 
locations along the lower-surface center line. These indicated transition locations 
taken from a relatively smooth thermocouple model were then compared with estimates 
determined from reference 12, which was a phase-change-paint study at a = 30°, to 
investigate the effect of simulated tile gaps on heating. The differences in transi- 
tion locations between the two studies were insignificant. The combined data were 
used to predict the flight-transition times represented by the vertical lines in 
figure 21 , with Reynolds numbers and angles of attack representative of the STS-1 
projected trajectory. 

At the forward fuselage location (x/c = 0.10), the two sets of flight data agree 
well with each other and with the estimates; there is no apparent indication of tran- 
sition at this location from the flight data at any time depicted during entry, with 
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the possible exception of the last point for STS-2 at approximately 1350 sec. The 
wind-tunnel data indicated probable transition at approximately 1390 sec. Transition 
is more evident at the downstream locations. For example, transition apparently 
occurred at approximately 1250 sec at x/c = 0.25 during both the STS-1 and STS-2 
entries. 

The wind-tunnel data indicated that transition would occur at approximately 
1360 sec. The laminar flight heating rates at this longitudinal location generally 
agree with the laminar estimates, and the turbulent flight heating rates agree with 
the turbulent estimates. This trend is generally true at other locations as well. 

The heating data at the 50-percent longitudinal location (x/c = 0.50) is noteworthy 
because of the differences in time at which transition appeared to occur during the 
two trajectories. Transition was apparently triggered a little after 1100 sec during 
the STS-1 entry, whereas it was delayed until almost 1300 sec during the STS-2 entry. 
The wind-tunnel data indicated probable transition at approximately 1270 sec. The 
flight heating rates prompted by the earlier transition are larger than the estimated 
laminar peak for the projected STS-1 entry, although they are well within the peak 
values of the aerodynamic-heating data book which are based on a nominal-design 
trajectory (14414.1). However, the extreme steepness of the estimated turbulent- 
heating curve indicates that if transition had occurred between 50 and 100 sec 
earlier, the heating rates would have exceeded the design-data-book curve as well. 

The transition times indicated by the wind-tunnel data for longitudinal locations 
greater than x/c “ 0.50 were much earlier than those experienced in flight. 


CONCLUDING REMARKS 

An assessment of the heating on the lower surface of the Shuttle orbiter was 
conducted at the Langley Research Center over the nominal-design entry trajectory for 
aerodynamic heating (14414.1 trajectory) and the proposed STS-1 (first flight of the 
space transportation system) trajectory. Perfect-gas laminar radiation-equilibrium 
heating rates and surface temperatures were calculated based on a combination of 
thermocouple and phase-change-paint wind-tunnel data. Averaged data from a large 
number of phase-change-paint tests at identical test conditions were successfully 
applied in this study to predict detailed quantitative heating results which were 
comparable in accuracy to thermocouple measurements. Therefore, the utilization of 
repeated phase-change-paint tests to form a suitable data base for extracting 
averaged quantities should not be ignored as a viable approach for effectively 
analyzing heating over complex configurations. 

The time histories of heating rates calculated for the nominal-design 14414.1 
trajectory along the center line of the orbiter lower surface and the 40- , 55-, and 
70-percent locations of maximum span were compared directly with curves contained in 
the aerodynamic-heating data book. Except for the 55-percent span location, the 
design data were approximately 40 percent higher during peak laminar heating than 
those determined in the present study. The design-data heating rates for the 
55-percent span location were less than 10 percent higher than the laminar-heating 
estimates of this study. 

Peak laminar heating rates during the proposed STS-1 trajectory tend to be of 
lesser magnitude but of somewhat longer duration than those estimated for the 
nominal-design 14414.1 trajectory. This result was not unexpected because for an 
insula tor- type thermal protection system (TPS) such as that of the orbiter, the 
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maximum temperature is the most significant parameter. The calculated laminar 
heating rates of this study compare favorably with flight data from the STS-1 and 
STS-2 trajectories. 

Turbulent heating rates were also calculated to represent a maximum heating 
boundary over the trajectories. Peak turbulent heating rates are estimated to occur 
200 to 400 sec later into the entry than the peak laminar rates and can exceed the 
design data by 100 percent or more. The calculated turbulent heating rates also 
compare favorably with STS-1 and STS-2 flight data over those portions of the 
trajectories where transition has obviously occurred. Typical turbulent heating 
rates occurring prior to 1200 sec into the entry of either trajectory generally 
exceed the corresponding laminar peaks; and prior to 1100 sec, the heating rates 
tend to exceed the design-data peaks. Therefore, the time at which boundary-layer 
transition takes place has a significant impact on the TPS performance of the 
orbiter. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
April 22, 1983 
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Figure 5.- Continued. 
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Figure 6.- Phase-change-paint heat-transfer-coefficient ratios ^/^^ef 
a = 35° and « 1 .79 x 10®. = 7.94; = 150°F; data taken 

from reference 4. 
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Figure 7.- Continued. 


32 
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T__ = 250° F; data taken from reference 5. 


Figure 7.- Continued 
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Figure 8.- Phase-change-paint heat- transfer-coefficient ratios 
U/U^g£ for a » 30° and R^ ^ » 7 x 10®. 
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Figure 8.- Continued. 
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Figure 8.- Concluded. 
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Figure 10.- Wind-tunnel experimental heat- transfer-coefficient ratios from references 3 and 4 
applied to nominal-design 14414.1 trajectory between 600 and 850 sec of entry. « 8; 

a « 35°; R « 2 X 10°. 
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Figure 11.- Wind-tunnel experimental heat-transfer-coefficient ratios from references 3 
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Figure 13.- Continued 
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Figure 16.- Pressure-coefficient comparisons from fuselage lower-surface center line. Tests were made 
in the AEDC Tunnel B and the Ames 3.5-Foot Hypersonic Wind Tunnel. 
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Figure 18,- Continued. 
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Figure 19.- Laminar and turbulent radiation-equilibrium heating rates for 
proposed STS-1 entry trajectory overlayed on equivalent set for 
nominal-design 14414.1 entry trajectory. 
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(b) x/c = 0.25. 
Figure 19.- Continued 
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Figure 19.- Continued 
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Figure 21.- Comparison of estimated and measured flight heating rates along lower-surface 

center line of Shuttle orbiter. 
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